Vector vortex solitons in nematic liquid crystals 
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We analyze the existence and stability of two-component vector solitons in nematic liquid crystals for which 
one of the components carries angular momentum and describes a vortex beam. We demonstrate that the 
nonlocal, nonlinear response can dramatically enhance the field coupling leading to the stabilization of the 
vortex beam when the amplitude of the second beam exceeds some threshold value. We develop a variational 
approach to describe this effect analytically. © 2008 Optical Society of America 
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Optical vortices are usually introduced as phase singu- 
larities in diffracting optical beams [1] and can be gener- 
ated in both linear and nonlinear media. The well known 
effect accompanying the propagation of such singular 
beams and vortex solitons in self-focusing, nonlinear me- 
dia is vortex breakup into several fundamental solitons 
via a symmetry-breaking azimuthal instability [2] . How- 
ever, recent numerical studies have revealed that spa- 
tially localized vortex solitons can be stabilized in highly 
nonlocal self-focusing nonlinear media [3-5]. This stabi- 
lization effect was later explained analytically [6] by em- 
ploying a modulation theory for the vortex parameters 
based on an averaged Lagrangian. 

Spatial optical vector solitons can form when several 
beams propagate together, interacting parametrically or 
via the effect of cross-phase modulation [7]. The sim- 
plest vector solitons are known as shape-preserving, self- 
localized solutions of coupled nonlinear evolution equa- 
tions [7]. A novel class of vector solitons in the form 
of two color spatial solitons in a highly nonlocal and 
anisotropic Kerr-like medium were predicted to exist in 
nematic liquid crystals [8-10]. The first experimental 
observations of anisotropic, nonlocal vector solitons in 
unbiased nematic liquid crystals were reported by Al- 
berucci et al. [10], who investigated the interaction be- 
tween two beams of different wavelengths and observed 
that two extraordinarily polarized beams of different 
wavelengths can nonlincarly couple, compensating for 
the beam walk-off, so creating a vector soliton. 

The main purpose of this Letter is twofold. Firstly, 
we introduce a novel class of vector solitons in nonlo- 
cal, nonlinear media, such as nematic liquid crystals and 
study their properties. These vector solitons appear as 
two color, self-trapped beams for which one of the com- 
ponents carries angular momentum and describes a vor- 
tex beam. Secondly, we demonstrate that the nonlocal, 
nonlinear response may dramatically enhance the field 



coupling, leading to the stabilization of the vortex for 
much weaker nonlocality when the amplitude of the sec- 
ond beam exceeds some threshold value. We develop a 
variational approach to describe this effect analytically. 

We consider the propagation of two light beams of dif- 
ferent wavelengths in a cell filled with a nematic liquid 
crystal. The light propagates in the z direction, with the 
(x, y) plane orthogonal to this. The electric fields of the 
light beams are assumed to be polarized in the x direc- 
tion. The system of coupled equations for the dimension- 
less complex field amplitudes u and v can be written in 
the form 
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where describes the change of the director angle from 
the pre-tilt state, which is related to the nonlinear correc- 
tion to the optical refractive index. In Eqs. |T]) the longi- 
tudinal (z) and transverse (x, y) coordinates arc normal- 
ized to the diffraction length and the beam width, re- 
spectively. The parameter v describes the degree of non- 
locality of the nonlinear response. When v — > 0, Eqs. ([lj) 
reduce to the Manakov vector nonlinear equations. The 
system (TTJ) conserves the energy flow P = Pi + P 2 = 
SS-™(W\ 2 + \v\ 2 )dxdy. 

We are interested in a special class of vector solitons 
for which one of the components carries angular momen- 
tum and the other component describes a spatially local- 
ized mode in the form of a spatial bright beam. Solutions 
of this type have been discussed earlier for nonlinear sys- 
tems with a local response [11,12] and they have been 
shown to be unstable in a large region of their existence 
domain [13]. For our system described by Eqs. |T]) such 
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solutions can be found in the form u = u>i(r)exp(i6iz) 
and v = W2(r)exp(i(j))exp(ib2z), where Wi(r) and W2(r) 
are real functions describing the beam envelopes, 
are real propagation constants and r = \/ x 1 + y 2 in 
the radial coordinate. The resulting system of equations 
obtained after substitution of these solution forms into 
Eqs. ([1]) is solved using a standard numerical relaxation 
method. Without loss of generality, we search for solu- 
tions with 62 < bi and set b\ = 3 to investigate the 
properties of vector vortex solitons by varying the prop- 
agation constant 62 and the nonlocality parameter v. 




Fig. 1. (a,b) Field distribution for (a) bright soliton and 
(b) vortex beam, (c) Nonlinear correction to the refrac- 
tive index, (d) Phase structure of the corresponding vor- 
tex beam shown in (b). Here v = 1, b\ = 3, and 62 = 1-8. 

Figure Q] presents an example of vector vortex soli- 
tons for Eqs. ([I]) for which one component has the shape 
of a bright soliton [Fig. [Ha)] and the other compo- 
nent carries angular moment, so forming a vortex soliton 
[Fig. [ljb,d)]. Due to the physical nature of the nonlocal 
response of the nematic liquid crystal, we notice that the 
refractive index change features a bell-shaped distribu- 
tion [Fig. [TJc)] , even though there is a singularity in the 
center of the vortex beam, this being crucial for the sta- 
bilization of vortex solitons. As shown in Fig. [U^a), for 
fixed propagation constant 61 and nonlocality parameter 
v the power of the vortex beam is a monotonically in- 
creasing function of the propagation constant 62, while 
the power of the bright soliton decreases monotonically. 
It is important to note that vector vortex solitons exist 
in a finite band of the propagation constant 62- At the 
lower band edge the vortex beam vanishes and one ob- 
tains a scalar bright soliton. However, at the upper band 
edge the beam with a bell shape vanishes, so that the 
vector vortex soliton transforms into a scalar vortex soli- 
ton. We find that the existence domain of vector vortex 
solitons shrinks with increasing nonlocality parameter v 
[Fig. Eft)]. 
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Fig. 2. (a) Power of bright (dashed) and vortex (solid) 
beams for v = 1. (b) Existence domain (gray) of vector 
solitons as a function of nonlocality parameter v (at b\ = 
3). 
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Fig. 3. Dependence of the critical value of the amplitude 
for a bright soliton on the nonlocality degree v. 



As one of the central results we find that a bright beam 
with a finite amplitude can stabilize an otherwise unsta- 
ble vortex beam. To address this issue we performed ex- 
tensive numerical simulations of Eqs. ([1]) using the beam 
propagation method. Firstly, we employed a stationary 
form of the vortex beam (in order to minimize radiation) 
as an input beam for the v component only (namely, 
there is no bright soliton as an input for the u com- 
ponent), noting that a vortex beam is unstable when 
propagating alone. Then we added a Gaussian pulse to 
the u component and studied the dynamics of a vortex 
beam propagating together with a Gaussian beam by 
varying the nonlocality parameter v. Our main results 
are presented in Fig. [3j from which one can see that for 
low nonlocality, a higher amplitude of the bright beam 
is required for stabilization of the vortex beam, while for 
high enough nonlocality, the vortex beam is observed to 
be stable, even when propagating alone. Figure [4] shows 
some illustrative examples. It is clearly seen that when 
the vortex beam propagates alone it becomes unstable 
and breaks up into two filaments [see Fig. BJb)] • How- 
ever, when we add a Gaussian beam with amplitude 0.9 
above the threshold (namely, a[ t = 0.9), then the vortex 
beam co-propagates with the Gaussian beam in a stable 
manner. Thus, we draw the conclusion that a nonlocal, 
nonlinear response can dramatically enhance the field 
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coupling, leading to the stabilization of the vortex soli- 
ton when the amplitude of the Gaussian beam exceeds 
some threshold value. 
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Fig. 4. (a,b) Input and output of the intensity for an 
unstable vortex beam propagating alone in the medium. 
(c,d) Stable propagation of the vortex beam coupled to 
the bright component. Here v = 4 and a J = 0.9. 

Let us now give a sketch of the modulation theory 
analysis that allows deeper insight into the effect of the 
vortex stabilization and which will explain in a quantita- 
tive manner the results shown in Fig. [3] Using the two- 
component equations re-written in a Lagrangian formu- 
lation [14], we employ the trial functions for the vortex, 
its soliton component and the director angle, 



u = a u sech(r/w u ) e l<T " : 

„ — r/w v id>+ia v z 

v = a v re ' e v 
8 = a u sech 2 (r//? u ). 
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These functions are substituted into the Lagrangian, 
which is then averaged by integrating in r and (j> over 
the plane. This procedure gives the averaged Lagrangian 
C = C u + C v + C uv , where C u is given in Ref. [15] 
and C v is given in Ref. [6]. The important interac- 
tion Lagrangian between the vortex and the soliton is 
Cuv = (a> 2 )(2V2^)" 1 a*v%. 

The results of Minzoni et al [6] show that a vortex is 
stabilized as its width decreases and its amplitude in- 
creases. Therefore, we just need to show that the width 
of the vortex decreases as the amplitude of the soliton in 
the other component increases. If A v = a v w v e^ 1 is the 
amplitude of the vortex, then from Ref. [16] 



p 2 A 2 "\ 

— -= W v + (a u w u ) W v - - = 0. 

oy v z 



(3) 



Using the vortex width determined by this expression in 
the stability threshold of Ref. [6] , we find that the vortex 
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Alwi < 14.4. 



(4) 



Combining this criterion with the amplitude- width rela- 
tion © , it can be seen that the vortex stabilizes for lower 
values of nonlocality parameter v as the amplitude a u of 
the bright soliton increases, which explains the stability 
results shown in Fig. [3] and Fig. [4] 

In conclusion, we have described theoretically a novel 
type of stable vector vortex soliton in nonlocal, nonlin- 
ear media, such as nematic liquid crystals. These solitons 
appear in the form of two color self-trapped beams for 
which one of the components carries angular momentum 
and is stabilized by the nonlocality enhanced beam cou- 
pling with the other spatially localized beam. We have 
studied the effect of this stabilization numerically and 
have also developed a variational approach to describe it 
analytically. 
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